Draft version February 2, 2008 

Preprint typeset using I^T^jX style cmulatcapj v. 6/22/04 



THE WARPED CIRCUMSTELLAR DISK OF HD100546 
Alice C. Quillen 

Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627 and 
Visitor, Research School of Astronomy and Astrophysics, Australian National University, Mount Stromlo Observatory, Cotter Road, 

Weston Creek, ACT 2611, Australia 
Draft version February 2, 2008 

ABSTRACT 

We propose that the two armed spiral features seen in visible Hubble Space Telescope images of 
scattered light in HD100546's circumstellar disk are caused by the illumination of a warped outer 
disk. A tilt of 6-15 degrees from the symmetry plane can cause the observed surface brightness 
variations providing the disk is very twisted (highly warped) at radii greater than 200 AU where the 
spiral features are seen. Dust lanes are due in part to shadowing in the equatorial plane from the 
inner disk within a radius of 100 AU. HD100546's outer disk, if viewed edge-on, would appear similar 
to that of Beta Pictorus. A disk initially misaligned with a planetary system, becomes warped due 
to precession induced by planetesimal bodies and planets. However, the twistedness of HD100546's 
disk cannot be explained by precession during the lifetime of the system induced by a single Jovian 
mass planet within the clearing at ~ 13 AU. One possible explanation for the corrugated disk is 
that precession was induced by massive of bodies embedded in the disk at larger radius. This would 
require approximately a Jupiter mass of bodies well outside the central clearing at 13 AU and within 
the location of the spiral features or at radii approximately between 50-200 AU. 
Subject headings: 



1. INTRODUCTION 

High angular resolution imaging of circumstellar disks 
have shown that these disks may not be coplanar. For 
example, the disk of Beta Pictorus exhibits a warp; 
the outer disk at radii greater than 50A U is tilted be- 
tween 3-5 degrees from the disk interior ijBurrows et alJ 
119951 iHean et al.l l20001. Models for the scattered opti- 
cal light of the bowed disk of AU Microscopii suggest 
tha t this disk is also warped, with a small tilt of about 
3° lIKrist et al .112001 . While warped disks are easiest 
to identify in nearly edge-on systems such as Beta Pic- 
torus and AU Microscopii, less inclined or nearly face- 
on disks may also not be coplanar. This leads us to 
search for features in less inclined disks that might also 
be explained with a war p ed disk model. As pro posed 
by iMouillet et aD l)1997|) : lAueereau et~aT] l)2001b|) . the 
twistedness of Beta Pictorus's warped disk could be a 
result of precession induced by unseen planets or plan- 
etesimals residing in the disk. Consequently by probing 
the geometry of the warp, we can constrain the proper- 
ties of the unseen planetary system. 

HD100546 is a nearby southern Herbig Be star 
(B9.5Ve; distance d = 103 ± 6pc with an ag e esti- 
mated t ~ lOMyrs l|van den Ancker et alJ 119971) which 
exhibits a dusty circumstellar disk revealed in scattered 
light from visible and near-in frared imaging ijGradv et alJ 
120011 lAugereau et al.ll2001a|) . We note that at the dis- 
tance of HD 100546, 1" corresponds approximately to 
100AU. Most noticeable in the images of this disk taken 
by the Hubble Space Telescope (HST) using the STIS 
camera are the two nearly symmetrical s piral features 
at a r a dius of about 250 AU from the star l|Cradv et alJ 
1200 1|) . iQuillen et all (|2005|) showed that the two arms 
were unlikely to be caused by spiral density waves ex- 
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cited by a planet embedded in the disk but could be 
density variati ons pulled out t idally by a nearby passing 
star. However, IQuillen et al.l 1)200 5) then ruled out this 
scenario for a number of reasons. To match the open- 
ness of the arms, the fiyby had to be recent; within a few 
thousand years. Not only is the probability for this sce- 
nario extremely low, but it implies that the perturbing 
object is near HD100546. A search for a perturbing ob- 
ject (star) mas sive enough to accou nt for the two arms 
was negative. IQuillen et alJ l)2005fl concluded that the 
spiral features were not density enhancements resulting 
from the tidal fo rce of a recent stella r clo se encounter. 
In con trast, both IQuillen et al.l l)2005[) and lArdila et alJ 
( 2005) have shown that the spiral arms in the outer part 
of HD141569A's disk are likely to haven been pulled out 
tidally by the nearby passage of it's binary companion, 
HD141569B,C. 

Warped disks can display complex morphology. 
Previous work modeli ng ga la ctic warped disk s (e.g., 
Quillen et all [ 19921 11993 IQuillen fc Bowerl [ 1999 1 : 
Nicholson et all Il992t iTubbsl fl980t IQuillen et alJ bo05t 
Steiman-Cameron et all I1992T) accounts for the optical 
dust lanes seen in Centaurus A and other galaxies and 
the parallelogram shaped infrared morphology seen in 
Spitzer images of Centaurus A. These previous studies 
have shown that the apparent morphology of a warped 
disk can be sen sitive to its optical properties (e.g., 
iKrist et alJl2005j) as well as its geometry. For emitting 
optically thin disks, folds in the disk can a ppear brighter 
to th e viewer than other regions (see IQuillen et alJ 
120051 and associated simulations). Likewise we expect 
that optically thin disks appear brighter in scattered 
light in regions where the disk tilts at high inclination 
with respect to the viewer. This is because the surface 
brightness of an optically thin medium of uniform 
thickness that scatters light from a distant source is 
dependent on its orientation angle with respect to the 



2 



viewer. 

The outer disk is illuminated by starlight that may pass 
through an inner absorbing disk. Consequently starlight 
illuminating the outer disk should be attenuated in the 
equatorial plane. A warped disk can be described as a 
series of rings. Each ring crosses the equatorial plane at 
two points. The shadow from the inner disk could cause 
two spiral dustlanes or shadows in the outer warped disk, 
resulting in the appearance of two spiral features. 

Here we consider the possibility that the disk of 
HD 100546 is warped, and that the warp is responsible 
for the observed spiral features. We first consider the 
case where the outer disk is isotropically illuminated by 
the central star. If the disk at large radii is warped and 
optically thick then we would expect the optical images 
would exhibit only one spiral arm. Only the high side 
of the disk would be illuminated whereas the opposite 
side would be in shadow. If the disk is optically thin 
and warped at radii greater than 2" or 200AU, then we 
expect to see two spiral features. Regions of high incli- 
nation with respect to the viewer would appear brighter 
and would be located on opposite sides of the star. 

We now consider the situation of a warped disk which 
is not isotropically illuminated. We expect that starlight 
illuminating the outer disk should be highly attenuated 
in the equatorial plane. At each radius an inclined disk 
crosses the equatorial plane at two opposing points. Con- 
sequently the shadow from the inner disk could cause two 
spiral dustlanes or shadows in the outer warped disk, re- 
sulting in the appearance of two spiral features. If the 
spiral features are due to an equatorial shadow then the 
outer disk could either be optically thin or thick. If the 
disk is optically thick then one would expect 2 spiral 
shadows and one bright spiral feature, wheres if the disk 
is optically thin then one would expect 2 spiral shadows 
with brighter spiral features lying between them. 

Unfortunately it is not straightforward to estimate the 
disk optical depth and width of HD100546's disk as a 
function of radius. Based on the near-I R scattered light 
surface brightness, IPantin et alJ l|2000|) estimated that 
the disk has a high normal optical de pth, r ~ 1, at 
radii gr eater than 1" or 100AU. However. I Aueereau et al.l 
(2001a) estimated a normal optical depth 40-50 times 
lower implying that the disk was optically thin at optical 
and near-IR wavelengths outside of ~ 80 AU. The ob- 
served infrared spectral energy distribution of HD100546 
implies that the luminosity of emission from dust is large, 
50% of the stellar photospheric luminosity. Most of the 
disk luminosity, 70%, is emitted in the mid-infrared an d 
so from a radius ~ 10 — 20 AU ijBouwman et alJ l20031. 
well within the location of the spiral features at ~ 250 
AU in scattered light. The estimated radius of the mid- 
infrared emission based on the spectral energy distribu- 
tion is approximately consistent with mid- i nfrare d high 
angular resolution observations (|Liu et alJ 12*003'). The 
infrared spectral energy distribution can be explained 
with a disk containing an inner, puffed up region ab- 
sorbing and re-radiating half of the stellar light, account- 
ing for the large mid-infrared flux. Because of its large 
covering a ngle the puffed inner disk edge shades the 
outer disk llBouw man et all2003UDullemond et af1l2001t 
iDullemond fc Dominikll2004[) . However, this model im- 
plies that the disk greater than 1" or 100AU which has 
been observed in scattered light is illuminated by the 



starlight that must pass first through the inner puffed 
up region of the disk. This situation makes it challeng- 
ing to estimate the optical depth of the disk from the 
scattered light at radii greater than 1". 

In this paper we search for warped disk models which 
can account for the spiral features seen in HD100546's 
disk. Warped disk models have some advantages over 
other models. Because planetary systems are nearly Ke- 
plerian, orbit orientations vary extremely slowly. Con- 
sequently a warp can be maintained for many rota- 
tion periods. This eliminates one of the problems of 
the previously proposed transient spiral structure model. 
Edge-on systems such as Beta Pictorus which are clearly 
warped, suggest that less highly inclined systems such 
as HD100546 could also be warped. In section 2 we de- 
scribe how we represent the geometry of a warped disk. 
We then explain how we synthesize model surface bright- 
ness images which can be compared to the observations. 
In section 3 we discuss dynamical or physical models for 
the disk warp based on the geometrical model which best 
matches the observed disk morphology. A summary and 
discussion follows. 

2. WARPED DISK MODELS 

2.1. Warp geometry 

We first discuss our notation for describing the orienta- 
tion of the planetary system with respect to the viewer. 
The orientation of a coplanar planetary system requires 
2 angles to describe; X> corresponding to the position 
angle (counter clockwise from North) of the axis of disk 
rotation on the sky, and an inclination angle, which 
describes the tilt of this axis with respect to the line of 
sight. If $ is zero, then the planetary system is viewed 
face-on. For a system that is not coplanar, \ and •& re- 
fer to the orientation of the rotation axis corresponding 
to the total angular momentum of the system's disk and 
planets. 

We describe the warped disk with respect to the ro- 
tation axis of the planetary system. A warped disk un- 
dergoing circular motion can be described as a series of 
rotating tilted rings, each with a different radius, r. The 
orientation of each ring is specified by two angles, a pre- 
cession angle, fi(r), similar to the longitude of the accent- 
ing node, and an inclination angle, i(r). These angles are 
given with respect to the rotation axis of the planetary 
system and the direction of the line of sight. We measure 
the angle f2 from a reference line in the system's ecliptic. 
As viewed on the sky, this line lies on the rotation axis 
of the system, but is projected onto the system's ecliptic 
plane. Our precession angle is the angle (measured at 
the star) between this line and the point at which the 
ring crosses the system's ecliptic. This is similar to the 
longitude of the ascending node which is measured with 
respect to the vernal equinox and the point at which the 
orbit crosses the ecliptic. 

2.2. Constructing a model scattered light image 

To produce an model image of the optical scattering 
light, all reflecting and absorbing regions along the line 
of sight at each position on the sky must be consid- 
ered. When the disk is optically thin, multiple scattering 
events and absorption can be neglected. In this case each 
photon originates from the star and is then reflected from 
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a single spot on the disk. At each position on the sky we 
can sum the reflected light at each location in the disk 
along the line of sight. We have restricted this modeling 
effort to an o ptically thin outer disk, consistent with the 
estimates by lAugereau et al.l l)2001aj) at radii outside of 
80AU. However we must keep in mind that future mod- 
eling efforts may need to consider outer disks with higher 
optical depth. 

We begin by randomly sampling x, y positions in the 
plane perpendicular to the disk rotation axis. At each 
position we compute a disk plane z coordinate based on 
our assumed function for the precession and inclination 
angles 51 (r) and i(r). To account for the disk thickness 
we add a vertical offset to z which was randomly chosen 
from a Gaussian distribution function. For the vertical 
structure of the disk we assume a normal distribution 



p(z) oc exp(— z 2 /2h 2 ) 



(1) 



where h is the standard deviation of the distribution. 
The FWHM of this distribution is 2Mh. We adopt a 
disk aspect ratio, h/r, that is independent of radius. 

The coordinates of the x, y, z position (in coordinates 
defined by the disk rotation axis) are then rotated using 
X and i? to account for the orientation of the system ro- 
tation axis with respect to the viewer. To produce each 
model image, 200,000 points in the disk are sampled. 
Scattered starlight from these points are summed along 
each line of sight to produce a surface brightness image 
on the sky. The brightness of each point in the disk de- 
pends on the assumed albedo times the normal optical 
depth of that portion of the disk multiplied by the flux 
from the star at that position. The scattering amplitude 
of the reflected light from each point in the 3D model 
was modified by the Henyey-Greenstein scattering phase 
function. The scattering asymmetric parameter, g, de- 
scribes the scattering anisotropy (g — corresponds to 
isotropic scattering, g = 1 to fully forward scattering). 
However, we do not vary g but instead adopt a fixed 
value of g = 0.5. We chose a representative value for 
g because of the uncertainty in the dust grain distribu- 
tion and composition and because of the broadness of the 
filter used for the HST observations. 

We take into account absorption of starlight from the 
inner disk by attenuating the starlight reaching the outer 
disk. We use an axisymmetric attenuation function that 
depends on the spherical coordinate, 9, where sin# = 
z/r. We assume that the opacity of the inner disk 



r d {9) =T O exp(-|0|/0 r ) 



(2) 



where To is the opacity in the equatorial plane, and 9 T 
describes an angular scale length. Approximately 50% of 
the total stellar luminosity is absorbed by the inner disk. 
An opaque torus blocking 50% of the light would cover 
angles -30° < 9 < 30°. Because the star is an B9.5V 
star its spectrum is quite blue. Consequently we expect 
T d 1 at broad band optical wavelengths for angles a few 
times smaller than the 9 ~ 30° required to block 50% of 
the total stellar luminosity. 

Because the disk surface brightness is observed to 
drop oc r~ 3 within 2.5 " of the star ijGradv et al J 120011: 
lAugereau et al.ll2001al) . the disk normal optical depth 
must depend on radius. Taking into account the r~ 2 
flux from the star, we match the observed radial drop 
with a —1 radial power in the normal optical depth 



ijAugereau et al.l2001aj) . The normal optical depth times 
the albedo is taken to be a power law function of radius, 
oc r -1 . We note that if there are large variations in 
the disk orientation then a different radial function may 
provide a better fit to the surface brightness profile of 
scattered light. 

2.3. The precession angle 

If the warp is due to the tidal force of a planet interior 
to the disk, then the precession rate of the longitude of 
the ascending node is approximately 
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as discussed by iMouillet et al.l l)1997|) . This approxima- 
tion is appropriate for particles at low inclinations and 
in nearly circular orbits. Here n is the mean motion of 
the disk at radius r (equivalent to the angular rotation 
rate for a circular orbit), D is the semi-major axis of the 
planet, M p is the mass of the planet, and M* is the mass 
of the central star. This precession rate is appropriate in 
the limit for D <C r and is independent of the disk incli- 
nation. This precession rate also neglects the self-gravity 
of the disk and so is only appropriate when the disk is 
low mass. 

After a time At, an initially flat disk will have a pre- 
cession angle 



(4) 



where fio is the precession angle at a reference radius, 
ro. When the precession is due solely to a single distant 
planet, the constant 
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and (3 = 3.5, corresponding to precession angle Q oc 
r~ 7 / 2 . However if there are multiple planets in the disk 
or the disk itself contains contains mass, one may con- 
sider more general laws or power laws with (3 < 3.5. In 
our numerical exploration we have explored variations 
in do, (3 and A to match the observed morphology of 
HD100546's disk. 

If a planet internal to the disk is initially taken out 
of the plane containing the disk then the tilt angle i 
(with respect to the plane containing the planet) de- 
cre ases with incre a sing r adius; the situation considered 
bv IMouillet et all <| 1997ft . However if the disk itself is 
tilted via tidal forces from an external stellar encounter 
then i would increase with increasing radius. Conse- 
quently we allow the disk tilt angle, i(r), to vary slowly 
with radius. At small radii where structure in the disk 
is difficult to resolve in the presence of scattered light 
from the star, models of the form given by Equation (0J 
predict extremely tight corrugations. We let i drop to 
zero at small radii, so that tight corrugations in the in- 
ner region were removed from the model images. This 
allowed us to keep a simple powerlaw form for the pre- 
cession angle. We allowed the inclination to vary with 
radius smoothly by using a spline function specified at 4 
different radii r = 0, 100, 200 and 400AU. The inclination 
at r = was set to 0°. 
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Common parameters for models are listed in Tabled 
Parameters varied for individual models are listed in Ta- 
ble |21 We first attempt to match the morphology of the 
disk with a geometric model and then discuss physical 
models which can account for the observed geometry. 

3. MODEL DISKS 

In Figure ^ we show a model warped disk in compar- 
ison to the STIS i mage HD100546 from Figure le by 
IGradv et alJ ((2001). The parameters used to describe 
this model are listed in Table 2] and as model #MA in 
Tabled The rings comprising the disk are projected onto 
the sky in Figure From a comparison between Figure 
Eland Figure G] we see that locations where the rings are 
in close proximity correspond to regions of higher sur- 
face brightness. These are regions where the disk slope 
or surface gradient (with respect to the line of sight) is 
high. A shadow from the inner disk is seen along the ring 
minor axes. 

Surface brightness profiles for Model #MA along the 
major and minor axes, at position angles 127 and 37°, are 
shown in Figure|3] We compare our model surface bright- 
ness p rofiles with the maj or axis profile shown in Figure 
5a by IGradv et a,U 1I200I of the STIS image. In both 
our model and the observed disk one spiral feature cor- 
responds to a bump in the major axis surface brightness 
profile on the northwestern side at about 300AU from 
the nucleus. The bump in the surface brightness profile 
is about 0.2 in the log above a smoothly dropping curve; 
approximately consistent with that seen in the observed 
profile. The increase in surface brightness corresponding 
to the opposite spiral feature is much less prominent on 
the southeastern side in both model and observed profile 
than that on the northwestern side. Dustlanes are not 
prominent in the major axis surface brightness profile (in 
both model and observed profile), however they are deep 
in the model minor axis profile. The dips in the major 
axis profile exhibited by our model in Figure Otj corre- 
sponding to dustlanes, are deeper than those observed. 
Equatorial shadowing in the model is more extreme than 
that observed. We discuss this problem in more detail 
below. The southern side of the disk is brighter than 
the north side because we allowed the scattering to be 
anisotropic and have taken the southern side to be nearer 
the viewer than the norther side. Some anisotropy in the 
scattering is consistent with the excess surface brightness 
in the near- infrared seen on the sou thern side at r ~ 3" 
reported bv lAugereau et all l)2001al) . 

We find that a warped model such as model #MA 
shown in Figures ^El can provide a good explanation 
for the spiral features observed in the disk of HD 100546. 
The warp causes apparent surface brightness variations 
along two tightly wound spiral features, with darker re- 
gions withi n them. Th e dark er regions were described 
as lanes bv IGradv et alJ l)2001l) . Two effects can account 
for regions of lower surface brightness: these regions are 
less inclined or more nearly perpendicular to the viewer 
(more nearly face-on) , or they lie in the equatorial plane 
and so can be in the shadow of the inner disk. Higher 
surface brightness regions are those that rise above the 
midplane shadow. They also correspond to regions that 
have a steeper surface gradient or slope in the disk due 
to corrugations in the disk. Low inclination evenly il- 
luminated warps can exhibit large spatial variations in 



the scattered light surface brightness, however this only 
occurs if the disk is highly corrugated or twisted and the 
disk is optically thin. 

In Figure 21 lower right-hand panel, we show a nearly 
edge-on model #ME, with parameters listed in Tables ^ 
and 121 This model is identical to model #MA, shown 
in Figures 1-3 except the system is more highly inclined. 
Model #ME can be compared to the edge-on disk of Beta 
Pictorus. Beta Pictorus is less twisted in its outer region 
compared to this model and does not ha ve a shadow in 
its mi dplane. The one planet model of Mouill et et alJ 
(1997) would lead us to expect that the disk is more 
highly twisted at smaller radii. In future we may con- 
sider the possibility that some of the unresolved struc- 
ture in Beta Pictorus's or AU Microscopii's inner disks 
at visible wavelengths might be explained with a warp 
that extends to smaller radii. We note that asymmetries 
between emission from one side of the disk compared to 
the opposite side would be introduced because the optical 
depth could be high at folds, and because of scattering 
asymmetry. 

3.1. Sensitivity to parameters 

In Figure 01 we show the effect of varying some of the 
parameters used to describe the disk. Model #MA has 
a disk inclination that increases with increasing radius. 
When the inclination is held constant, as is true in Model 
#MC, shown on the top-right in FigureQJ the spiral arms 
do not extend as far to the south-west and north-east 
as is seen in the scattered light image. We find that 
the angular extent of each spiral feature is smaller than 
that observed if i is held fixed. If the inclination with 
respect to the system axis increases with radius, then 
the spiral features extend over a larger range of angles, 
as shown in Model #MB on the upper left in Figure 21 
Model #MB has an even higher outer inclination than 
Model #MA. Disks with larger tilts (larger i) tend to 
produce higher surface brightness variations in the spiral 
arms. However if i is increased past 10 degrees then the 
disk can be folded with respect to the viewer and this 
reduces the contrast of the spiral features. Model #MA 
and #MB have regions where multiple folds of the disk 
are encountered along the line of sight (see Figure^) and 
only the fold edges are regions of high surface brightness. 

It is interesting to note that the apparent ends of the 
two spiral arms are about 180° apart and oriented nearly 
along the major axis of the disk. Such a situation arises 
naturally from the warped models. In contrast, spiral 
arms that are due to spiral density waves are not ex- 
pected to end at locations 180° apart. The folds of the 
warped disk cause high surface gradient regions (corre- 
sponding to higher surface brightness regions) on either 
side of nucleus. However along the major axis these folds 
are oriented along the line of sight. As a result, spiral 
features in the model scattered light images tend to end 
along the disk major axis. Higher resolution models show 
a nested series of self-similar spiral features inside the 
outermost ones. This follows since we have adopted a 
power law form for the precession angle. 

Similar morphology to that observed is seen when the 
inner disk has a higher equatorial opacity tq but a shorter 
opacity angular scale length 8 T . Model #MD shows such 
a model. This model has a slightly larger outer disk scale 
height which has the effect of smoothing the model sur- 
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face brightness image. Had we left the disk scale height 
similar to that of model #MA this model would have 
had extreme contrast in its dust lanes. We find that 
thinner disks have more sharply defined spiral features 
and deeper dustlanes. When the angular scale length 
of the equatorial shadow is shorter, (6 T is smaller), the 
dustlane is more sharply defined. Higher equatorial inner 
disk opacity (large t ) causes deeper shadows. 

As we commented above, we suspect that Model #MA 
has deeper dustlanes than observed. The equatorial 
shadow for Model #MA has an opacity of 1 (at visi- 
ble wavelengths) at equatorial angles 8 = ±13°. This 
opacity is sufficiently high that it is approximately con- 
sistent with the absorption of ~ 50% of the stellar light 
from the inner disk. However this high opacity and large 
angular scale length (in the inner disk opacity) implies 
that the outer disk is illuminated by visible light that has 
been significantly attenuated by the inner disk; much of 
the outer disk has a tilt lower than 13°. Our model com- 
putes the surface brightness at one wavelength only cor- 
responding to a central optical wavelength for the broad 
STIS image. It may be possible to improve the model 
image by integrating and summing images at different 
wavelengths. It is also possible that wings of the stel- 
lar point spread function have smoothed the appearance 
of the outer disk, reducing the actual surface brightness 
contrast. The depth of the dustlanes in our model images 
can be decreased by increasing the angular scale length of 
the shadow, 9 T . However then the model midplane opac- 
ity must be reduced to reproduce the observed morphol- 
ogy. A midplane opacity below 1 would be unrealistic as 
we expect the midplane to be dense and optically thick 
at visible wavelengths. To improve the model we sus- 
pect that we would require a more complex function for 
the shadow than that given by Equation . This func- 
tion would necessarily be described by a larger number 
of parameters. A more complex model is difficult to con- 
strained with the STIS image but could be constrained 
by future high quality and multi- wavelength images. 

The twistedness of the disk is set by Aq. Higher \Ao\ 
corresponds to more highly wound spiral features. How- 
ever, variations in Aq can also change the radius of the 
spiral features if the radius tq is not simultaneously ad- 
justed. Because the spiral features are tightly wound we 
could not place constraints on the parameter which 
sets the dependence of a on radius, consequently we 
set = 3.5, consis tent with the model explored by 
Momllet et, alJ l)1997j) . The contrast between the surface 
brightness on the near side compared to that on more dis- 
tant side is larger if the scattering is more anisotropic; 
(the scattering asymmetry parameter g is larger). 

Based on our exploration of models with different pa- 
rameters we have found the following: Only models with 
highly twisted disks have tightly wound spiral features 
similar to those observed. Only models with inclination 
increasing with radius have spiral features that extended 
over a sufficiently large range of azimuthal positions. 
Only relatively thin disks, h/r < 0.15, have sufficiently 
sharp or fine features to be similar to the observations. 
We find a degeneracy between the disk inclination, the 
angular form of attenuation from the inner disk, and the 
disk thickness, because these parameters all affect the 
contrast or amplitude of the spiral features. We did not 
find models with morphology similar to that observed 



with disks higher than i ~ 15° and lower than i ~ 6°. 
Low inclination disks failed to exhibit sufficiently high 
surface brightness variations and high inclination disks 
exhibited multiple folds along the line of sight, reduc- 
ing the surface brightness variations. Models with small 
shadow angular scale lengths [6 T < 5°) have dustlanes 
that are excessively deep, suggesting that the upper lay- 
ers of the inner disk have an opacity distribution (as 
a function of 9) with a moderately large angular scale 
length. 

3.2. Mass constraints 

By matching the observed morphology of the disk with 
our model, we can estimate the extent that the disk 
is twisted. This is described by parameters Aq, and 
(see Equation We consider here the hypothesis that 
the twist was caused by an unseen inn er planet that 
is misaligned with the oute r disk (e.g., iMouillet et alJ 
ll997tlAugereau et al.l200 1b). We remind the reader that 
Equation J^J refers to to a tilted disk with negligible mass 
which is perturbed by an inner planet. The parameter 
Aq depends on the time since the disk was initially tilted, 
At, and the mass and semi-major axis of the hypotheti- 
cal inner planet causing the precession. We can assume 
that At is less than the age of the star or At < 10 7 years. 
This allows us to place a lower limit on the mass of the 
planet times the square of the planet semi-major axis. 
Using Equation |(2J| and replacing At with the age of the 
system, t age , we find 

Computing these quantities for our value of A (in radi- 
ans) and reference radius tq (listed in Table |2J), we find 

D 2 M P > 2.9 x \tfAU 2 Mj (7) 

where Mj is the mass of Jupiter. However, this con- 
straint is impossible to satisfy for a single Jovian mass 
at small radius. This constraint cannot be satisfied for 
a Jovi an mass plan e t with in the lit edge of the disk at 
13 AU lGradv et alJ <|2005[) . A brown dwarf sized object 
at D ~ 200AU would open a gap in the disk that would 
be observable in the images. The simplest scenario of an 
initially misaligned low mass disk and and a single inner 
planet fails to account for the twistedness of the disk. 

One way to account for the highly twisted disk would 
be if the disk itself contained significant mass. For exam- 
ple, a Jovian mass of planetesimals between 100-200 AU 
but misaligned with the outer disk could account for the 
twisted nature of the disk. If this mass is not confined 
to one body but extended then the radial power of the 
precession angle would be reduced. This might account 
for the somewhat better match of our model #MB with 
= 2.5 instead of 3.5. HD100546's disk appears to be 
more twisted than Beta Pictorus's disk, however Beta 
Pictorus is 10-20 times older so the extent of the twist is 
less of a constraint on the planetary and disk system. 

3.3. Leading or trailing arms 

It may in future become possible to measure the sense 
of rotation of HD 1 00546 's disk. The disk rotation axis 
could be similar to the star's rotation axis. Unfortunately 
the rotation axis of the star HD100546 is not known. 
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iClarke et alJ 1)1999(1 attempted to measure this axis from 
polarization measurements, however these were affected 
by scattering from the surrounding dust. In spite of this, 
we now discuss the sense of the warp in comparison to the 
direction disk rotation. Future measurements of cither 
the disk or stellar rotation could be used to support or 
refute dynamical explanations for the spiral morphology. 

Equation shows that the inner disk should pro- 
cesses faster than outer disk and in the retrograde sense. 
Consequently we expect the disk to should twist in the 
direction of rotation as the radius increases. If we dif- 
ferentiate Equation we find that dSl/dr is positive 
and so Q increases in the direction of rotation. The spi- 
ral features would be leading instead of trailing. If the 
twist is a result of precession induced by bodies in the 
disk interior, we would predict that the disk is rotating 
clockwise on the sky. We expect the northwestern side to 
be blueshifted from the system center of mass, and the 
southeastern side to be redshiftcd. The southwestern side 
is probably closer to us, based on t he nebulosity on the 
southern side seen in the STIS image lGradv et all ( 2001) 
and because the excess seen in the scatter ed near-infrared 
light a t r ~ 3" on the southern side l)Augereau et al.l 
200l3. 

4. SUMMARY AND DISCUSSION 

In this paper we have presented a new model for the 
spiral features seen in visible scattered light images of 
HD100546's disk. We reproduce the two-armed spiral 
features with a highly twisted warped disk model. The 
disk inclination with respect to the systems ecliptic is 
low, 6 to 15° degrees and increasing with radius. The 
disk is fairly thin with a ratio of FWHM/r of ~ 0.15. 
Surface brightness variations are due to a high surface 
gradient with respect to the line of sight in the folds of 
the disk. Dark lanes increase the contrast of the spiral 
features and are caused by equatorial shadowing from 
an inner unresolved disk. As a nearly Keplerian sys- 
tem can maintain a warp for long periods of time (secu- 
lar timescales rather than rotational periods), this model 
has some advantages over transient spiral models. The 
observed spiral features end at positions approximately 
located on the disk major axis (as seen on the sky). This 
is a feature common to moderate inclination warped disk 
models that would not be exhibited by a spiral density 
wave model for the spiral features. 

However, the morphology requires the disk to be so 
twisted that the tidal force from a Jovian mass object 
in the inner disk clearing cannot have induced the twist 
during the lifetime of the star. One possibility is that 
the disk warp could have been induced by a significant 
mass (Jovian mass) of objects that are inclined with re- 
spect to the outer disk, well outside the inner clearing, 
at intermediate radii between 100-200 AU. 

We now discuss proposed mechanisms for accounting 
for the tilt between the disk and objects in the inner disk. 
A stellar encounter could tidally induce a tilt in the outer 
disk. In this case we would expect the disk tilt or incli- 
nation would increase w ith radius (e.g., as calculated by 
iKobavashi k, Idal l2001). Our models better match the 
angular extension of the spiral feature when the tilt in- 
creases with radius, suggesting that this may be the case 
for HD100546's disk. Such a stellar encounter, while un- 
likely in the field, could have been probable in a denser 



stellar environment such as the star's birth cluster. An 
alternative scenario is that resonances between multiple 
pla nets cause an inclination incr ease in one of the plan- 
ets l|Thommes &: Lissauerl 12003'!. In this case, the disk 
would have tilt or in clination that decreases with radius 
l|Mouillet et alJll997|) . Also this model would likely re- 
quire vertical resonances for planets at large semi-major 
axes to account for the highly twisted disk. The differ- 
ent dynamical scenarios predict different radial forms for 
the angles used to describe the warp, i and f2. Conse- 
quently better observations and accompanied modeling 
of the HD100546's disk and other circumstellar disks will 
produce better constraints on the disk tilt and precession 
angle as a function of radius. These in turn will allow 
better tests of the dynamical models for the warp for- 
mation and evolution. The dynamical models have the 
capability of providing unique constraints on the mass 
distribution in the outer disk. 

In this paper we have used a very simple model to take 
into account equatorial shadowing due to the opacity of 
the inner disk. Unfortunately we find a degeneracy in our 
models between the disk tilt, the angular form of attenu- 
ation from the inner disk, and the disk thickness. This is 
because these parameters all affect the surface brightness 
contrast or amplitude of the spiral features. Shadowing 
is likely to be a strong function of color, consequently 
multi-color imaging may be able to probe the structure 
of the inner unresolved disk as well as better constrain 
the geometry and structure of the outer disk. Warped 
disks seen at visible wavelengths should exhibit asymme- 
tries due to optical depth variations, and asymmetries 
due to forward scattering. Both of these affects should 
be more severe at bluer wavelengths, providing a pos- 
sible way to discriminate between geometrical models. 
We note that some asymmetries and perceived clumps 
in visible light i n edge-on warped disks (e.g., as seen in 
AU Microscopii: lLiiJl2004[l . might in future be explained 
with a warped disk rather than a model containing ec- 
centric rings. However asymmetries observed in thermal 
or mid-i nfrared emission, suc h as is observed in Beta Pic- 
torus bv lTelesco et al.l l|2005D . could not be explained via 
a warped disk model as at these wavelengths we expect 
the disks to be more nearly optically thin the emission 
should be nearly isotr opic. 

The observations of lEisner et al.l l|2004|) suggest that 
HeBe stars are not highly warped. However, moderate 
inclination warps, such as found here would not violate 
the comparison of submillimeter and near-infrared po- 
sition angles and ellipticities. It is possible that highly 
tilted warped disks may provide alternate explanations 
for excess far-infrared emission in some systems as such 
disks cover a larger solid angle and so can absorb more 
stellar flux at large radii than a coplanar disk of the same 
thickness. 
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Fig. 1. — The left-hand panel shows the log surface brightness of model warped disk, # MA with parameters are given in Tables and 
|21 The middle panel shows log of the surface brightness of the HST image of HD100546 from Gradv et al. ( 2001). Our model shows that a 
warped disk can exhibit surface brightness variations in a two armed spiral pattern. The surface brightness variations are due to the slope 
or gradient changes in the disk surface with respect to the line of sight and an equatorial shadow from an inner disk. 
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Fig. 2. — Projected circular disk rings at evenly spaced radii corresponding to the disk shown in figure^and described by model #MA. 
Regions where projected rings are in close proximity correspond to regions of bright surface brightness (see figure 0. This is because these 
regions have large slopes or surface gradients with respect to the line of sight. The equatorial shadow lies along the minor axes of the rings. 
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Fig. 3. — Major and minor axis surface brightness profiles for model #MA, that is also shown in Figures and [5] Dips in the surface 
brightness profile are seen where there are dark dustlanes due to equatorial shadowing by the inner disk. Along the major axis at about 
300 AU, an increase in the surface brightness is seen at a fold in the disk. A hole at r = 13AU has been placed near the star, primarily to 
limit the flux range covered by the plot. For the major axis profile, positive r refers to southeastern side of the disk. For the minor axis 
profile positive r refers to the near or southwestern side of the disk. The y-axis shows log 10 of the surface brightness normalized to the 
peak value. 
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Fig. 4. — The effect of varying warped disk parameters. These models are similar to that (Model #MA) shown in FiguresQ|2] andl^l 
Varied parameters are listed in Table |21 On the top left is Model #MB which has a higher inclination at large radii than model #MA. 
The top right panel shows Model #MC which has constant inclination with radius. The bottom left panel shows Model #MD which has 
a higher equatorial opacity for inner disk (larger tq) than Model #MA but a lower inner disk opacity angular scale height (9 T ). On the 
lower right we show a nearly edge-on disk, model #ME. to illustrate the extent of the twist, and for comparison to edge-on systems such 
as Beta Pictorus. 
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TABLE 1 

Common Parameters describing the warp models 



Parameter Value 

X 40° 

A -260° 

r 250 AU 

n 230° 

P 3.5 



Note. — The models shown in Figures HI |21 151 and Bl have these parameters in common. The position angle of the system rotation axis 
on the sky is denoted by \. The precession angle at the reference radius ro is Qrj- The reference radius ro is given in AU (100AU ~ 1" on 
the sky for HD100546). The parameters Ao and /3 describe the sensitivity of the precession angle Q with radius or the twistedness of the 
warp; see Equation Hi . The angles x> ^0 and Aq are given in degrees. 



TABLE 2 

Parameters describing the individual models 



#MA #MB #MC #MD #ME 



■& 52° 52° 52° 52° 89° 

i 3,6,15° 3,5,20° 6,6,6° 3,6,15° 3,6,15' 

h/r 0.05 0.05 0.05 0.08 0.05 

t 3.0 3.0 3.0 5.0 3.0 

9 T 12.0° 12.0° 12.0° 5.0° 12° 



Note. — The inclination of the system rotation axis with respect to the viewer is i9. The inclination of the warped disk with respect 
to the system rotation axis is i. The angles and i are given in degrees. The disk aspect ratio, assumed to be constant with radius is 
denoted by h/r. The standard deviation of the vertical density distribution is h and the FWHM is 2.35/i. Inclination angles are given at 
three radii, r = 100, 200 and 400AU, and we set i = at r = 0. A spline function was fit between these values. The opacity to describes 
the opacity of the inner disk in the equatorial plane. The parameter 9 T describes the angular scale height of the inner disk's opacity in 
degrees. Surface brightness profiles, projected rings and images for model #MA are shown in Figures jz] |3] Images are shown in Figure 
Hfor Models #MB-#ME. 



